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We have developed a numerical procedure to clarify the critical behavior near a quantum phase 
transition by analyzing a multi-point correlation function characterizing the ground state. This work 
presents a successful application of this procedure to the string order parameter of the 5=1 XX Z 
chain with uniaxial single-ion anisotropy. The finite-size string correlation function is estimated 
by the density matrix renormalization group method. We focus on the gradient of the inversed- 
system-size dependence of the correlation function on a logarithmic plot. This quantity shows that 
the finite-size scaling sensitively changes at the critical point. The behavior of the gradient with 
increasing system size is divergent, stable at a finite value, or rapidly decreases to zero when the 
system is in the disordered phase, at the critical point, or in the ordered phase, respectively. The 
analysis of the finite-size string correlation functions allows precise determination of the boundary of 
the Haldane phase and estimation of the critical exponent of the correlation length. Our estimates 
of the transition point and the critical exponents, which are determined only by the ground-state 
quantities, are consistent with results obtained from the analysis of the energy-level structure. Our 
analysis requires only the correlation functions of several finite sizes under the same condition as a 
candidate for the long-range order. The quantity is treated in the same manner irrespective of the 
kind of elements which destroy the order concerned. This work will assist in the development of a 
method to directly observe quantum phase transitions. 

PACS numbers: 75.10.Pq, 75.40.Cx, 75.40.Mg 



I. INTRODUCTION 

Quantum phase transitions originating from quantum 
fluctuations have been extensively studied as a hot issue 
in condensed-matter physics. Several interesting charac- 
teristics of the transitions appear in the low-energy be- 
havior of the systems. Two types of approach can capture 
the phase transitions and critical phenomena precisely 
when the transition is continuous. One is analyzing the 
energy-level structure. The other involves considering 
the ground-state behavior. 

In the former approach, a standard method is to an- 
alyze the structure of the energy levels of finite-size sys- 
tems based on the finite-size scaling (FSS) assumption. 
For example, the scaled energy gap [l| is often used to 
estimate the boundary of the gapped phase as the tran- 
sition point. This is called phenomenological renormal- 
ization group (PRG) analysis. However, it is difficult to 
estimate the transition point when a logarithmic correc- 
tion appears in the dependence of the energy difference. 
A typical example is the Berezinskii-Kosterlitz-Thouless 
(BKT) Transition 0. To resolve this difficulty, the level- 
spectroscopy method has been developed [|[ and precise 
determinations of phase transitions have been success- 
fully made for various transitions in many models. Unfor- 
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tunately, this analysis is complicated in that appropriate 
adjustments of the procedure are required according to 
the type of phase transition, which must also be known 
in advance. 

In the latter approach, on the other hand, quanti- 
ties that characterize the ground state are carefully ob- 
served. One of these quantities is the multi-point corre- 
lation function. The long-range behavior of correlation 
functions shows whether the system exhibits long-range 
order. If a correlation function survives to be nonzero in 
the long-range limit, it is an appropriate order parame- 
ter. However, it is not easy to capture a phase transition 
using this strategy because reliable and precise data on 
correlation functions are necessary for large systems. The 
system sizes that are treated in numerical-diagonalization 
calculations are insufficient. For this reason, the latter 
approach has been employed in only a few studies. There- 
fore, no systematic procedure for analyzing ground-state 
quantities to capture quantum phase transitions has been 
established to date. 

In this paper, we develop a procedure to determine 
the transition point and critical exponents by analyzing 
correlation functions based only on the scaling assump- 
tion. A feature of this approach is that only the com- 
mon quantities under the same condition are treated ir- 
respective of the type of phase transition. We call the 
procedure ground-state phenomenological renormaliza- 
tion group (GSPRG) analysis. To confirm its validity 
and usefulness in detecting phase transitions, we apply it 
to a nontrivial ground state in the AF S = 1 XX Z chain 
with uniaxial single-ion anisotropy by the density matrix 
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renormalization group (DMRG) method 

In the isotropic case of this system, there exists a 
nonzero energy gap between the unique ground state and 
the first excited state, called the Haldane gap [10]. It is 
known that when anisotropy is introduced, of the single- 
ion type or of the XXZ-type exchange interaction, the 
Haldane gap decreases and finally closes. The region 
where the nonzero Haldane gap exists is called the Hal- 
dane phase. The phase diagram of the AF 5 = 1 chain 
with anisotropy of the two types, including the Haldane 
phase, has been extensively studied by analyzing the 
energy-level structure, assisted by the level-spectroscopy 
method @,i,[T3, El. 

It is well known that in many AF spin systems, the 
standard spin-spin correlation function gives so-called 
Neel order. In the ground state in the Haldane phase, 
however, the spin-spin correlation function decays ex- 
ponentially with a finite correlation length and the Neel 
order no longer exists. In this sense, the Haldane phase is 
a disordered phase. However, the string order parameter 
is known to characterize the ground state in the Haldane 
phase, in which the longitudinal and transverse string or- 
der parameters are nonzero [l2l [l3| . From the viewpoint 
of the string order, it is possible to treat the Haldane 
phase as an ordered phase and capture the phase tran- 
sition at the boundary. This approach has been applied 
to numerical-diagonalization data of the Haldane phase 
in 5 = 1 systems [a [H], [HI- Unfortunately, only very 
small system sizes were available and hence it was quite 
difficult to capture precisely the critical behavior of the 
string order near the transition point. 

In this situation, we can obtain numerical data of the 
string order of this model for much larger system sizes 
by using DMRG for GSPRG analysis. Consequently, it is 
possible to examine the phase transition at the boundary 
of the Haldane phase from the viewpoint of the string or- 
der. We compare our results with those from the analysis 
of the energy-level structure. This comparison provides 
a systematic and consistent understanding of the phase 
transition. 

This paper is organized as follows. The model Hamil- 
tonian and order parameters are defined in section [Til 
The analysis procedure which we have developed is in- 
troduced in section Hill The numerical results and dis- 
cussions are given in section IIVI Section [V] consists of a 
summary of this work and some remarks. 



II. MODEL HAMILTONIAN AND ORDER 
PARAMETERS 

We consider the following Hamiltonian: 

JV-l N 

n=J2 \J(s*s* +1 +svg! + J+j z sist + ^if£s?, (1) 



i=l 



1=1 



where N is the system size, Sfia = x,y,z) are spin-1 
operators, J and J z represent the XXZ-type anisotropic 



exchange interaction, and the parameter D represents 
uniaxial single-ion anisotropy. In this paper, energies are 
measured in units of J and hence we take J = 1 here- 
after. The boundary condition of the system is open. An 
antifcrromagnetic chain is usually characterized by the 
Neel order parameter defined as 



(2) 



where 0^ el («, j) is the Neel correlation function given by 



os«(*\i) = (-iy- j (s?s?). 



(3) 



Here (B) represents the ground-state expectation value 
of an arbitrary operator B. In the Haldane phase, the 
Neel order parameter vanishes. However, the string or- 
der introduced by den Nijs and Rommelse [12| appears 
instead. The string order parameter is given by 



0^= lim 0: tI (i,j) 

\i— j\— >oo 



(4) 



where the string correlation function 0" tr (i, j) is given by 



5? exp 
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Kennedy and Tasaki [l3| extensively studied the string 
order and applied a nonlocal unitary transformation to 
the Hamiltonian ([TJ. Thereby they obtained 



^str 



Og m (n) for a 



where Of CI 



lim i.. 



(6) 



and TL is obtained 



\i— j\— tea \ J i J j )i 

by applying the transformation to the original Hamilto- 
nian H.. In the transformed system, Z2 x Z2 symmetry 
emerges. Breaking of this symmetry is described by the 
behavior of the order parameters Of erlo (H). When the 
system is in the Haldane phase, this Z2 x Z2 symmetry 
is fully broken. When the chain is in a phase other than 
the Haldane phase, the full symmetry or a part of the 
symmetry survives [1, 0, [l6| • 



III. ANALYSIS PROCEDURE 

In this section we introduce our analysis procedure. 
The procedure consists of three steps. The first is to cal- 
culate the longitudinal and transverse order parameters 
by the DMRG method based on the finite-size algorithm 
[j, Q ■ Knowing the behavior of the order parameters en- 
ables us to observe the boundary of the Haldane phase 
briefly and to confirm that the transition at the bound- 
ary is continuous. If we apply the FSS analysis of the 
order parameters (l7j . we can roughly estimate the crit- 
ical point and exponents. However, there remain finite- 
size effects in these estimates. In order to eliminate the 
effects, we carry out a finite-size extrapolation derived 
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from the FSS formula. As the second step, we introduce 
a finite-size quantity which we calculate from the string 
order parameters. This quantity reaches the critical ex- 
ponent just at the phase boundary when we take the limit 
N — > oo. By examining the behavior of this quantity, it 
is possible to obtain the phase boundary and the critical 
exponent of the order parameter at the transition point. 
At the third step, we estimate the critical exponent of 
the correlation length near the transition point by ex- 
trapolating a finite-size quantity at the above transition 
point. The present procedure gives consistent and precise 
estimates for the critical point and exponents. 



A. Calculation of order parameters 

To calculate the order parameters, we use the finite 
DMRG algorithm with the acceleration algorithm intro- 
duced by White 0, 1, El. A correlation function, such 
as Eq. ([3]) or Eq. ([5]), is calculated as follows. First, we 
obtain a variational wave function of the ground state 
represented by the matrix product state (MPS): 

\S>) -Vrff/M^f/P]^ ..[A72-1W/2-1 

{«} M 

A s N/2i s N/2 + l 
-^■"JV/2-1i q JV/2 + 2 " Q ]V/2 + 2> Q JV/2 + 3 

V [N-2],s N ^ 2 y[N-l],s N ^ 1 \ I g s \ 

v aN-2,a N -! v a N -i,s N ^1*1*2 °N/ j 

({s} = si, s 2 , • • • , s N , {a} = a 2 ,a 3 , • • ■ , ajv-i), 

(7) 

where U and V are matrices satisfying 

U®*U® = 1, (» = 2,..-,JV/2-l), (8a) 

yWyVH = i j (j = N/2 + 2, • • ■ ,N —1), (8b) 

respectively, and A is the ground state of the renormal- 
ized Hamiltonian which is obtained by applying the den- 
sity matrix renormalization transformation to the Hamil- 
tonian ((T|). A truncation error e due to the cut off m-p 
is estimated by 1 — Tr (A* A), where niF is the number 
of states preserved in the DMRG iterations Q. Using 
this wave function, we estimate the expectation value 
°stAhj) as follows: 

0« tr (i, j) = - <$| S?5f +1 • • • |$) 
= -Tr(EW 1 EW..-Et 1] E®E [ £ 1] ••• 

E [N/2-l\ E [A\ _ E l - N/2+2] 

^N/2-l ^ N/2^N/2 + l ^N/2+2 

e% 1] e [ £le [ 1 +1] ■ ■ ■ e[ n - 2] e[ n - 1] ) . (9) 

Here, 1 is the identity matrix in spin-f space and S a is 
defined as S a — exp(iirS a ). Each £ is a matrix for an 



arbitrary local operator O defined by 

E o.o= E HoKX*2|o|4) 

S\ ,Sj ,S2 ,s' 2 

(10a) 

E o = E (*i ^ uti] ' Si ® ( u[iiA )* - 

(3 < i < N/2-1), 

(10b) 

E o ] o= E {s k \0\ S ' k ){ Sk+l \0\s' k+1 ) 

s k ,s' k ,s k+1 ,s' k+1 

(fc = N/2), 

(10c) 

s o = E ( s i\ 14) FL7l,S ' ® (V b1 ^)* , 

(Ar/2 + 2<j <iV-2), 

(lOd) 

E l o.o 1] = E (^-i|0|*WW|0|^) 

S JV - 1 , 8 ' N _ x , S TV i S 

(10e) 

When calculating a two-point correlation function such 
as Eq. ([3]), we replace all the operators S a with 1. 
We also choose an appropriate correlation function as 
the longest-ranged component from among 0^ tr (i,j) or 
O^eel f° r a fixed N under the open boundary con- 
dition. Three desirable conditions should be met, as fol- 
lows. 

1 . The measurement points i and j are as far as pos- 
sible from the edges. 

2. The correlation distance of \j —i\ is as long as pos- 
sible. 

3. The distance \j — i\ should increase in proportion 
to the system size N. 

In order to satisfy the above conditions, we take i = 
N/3 + 1 and j = 2N/3. Thus, we consider the order 
parameter 

O a {N, D, J z ) = O a {N/3 + 1, 2AT/3, D, J z ), (11) 

where a = x,z, O a represents O^eei or ^str- Note here 
that N/6 should be an integer. We emphasize again that 
Eq. (1111) i s useful in the case of the open boundary con- 
dition [19|. 

The error of the order parameters is estimated as fol- 
lows. 
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1. (O a ) mi is calculated in the case of mi = 0.95 x mp. 

2. (0) m _ m and the truncation error e are calculated 
in the case of m 2 = mj . 

3. The error estimation of the order parameters 
(O a )"™ r is defined by the following formula: 

(O a )Z7 = max [| (° a ) m2 - (O a ) mi |,|e (O a ) m2 1] . 

(12) 



In this paper, all numerical data have a truncation error 
e < 10- 7 . 



B. Finite-size scaling analysis 

In the general theory of phase transitions, the treat- 
ment of physical quantities depends on whether the tran- 
sition is continuous or discontinuous. If the transition is 
continuous, the critical behavior of bulk quantities is ex- 
tracted through the FSS analysis of finite-size quantities 
B- As we observe later, the DMRG data of correla- 
tion functions of the finite-size systems are continuous 
near the boundary of the Haldane phase. The exact- 
diagonalization data of the string order parameters are 
also continuous. Thus, it is possible to perform FSS anal- 
ysis of our DMRG data of the string order parameter. 

The present model includes two control parameters, D 
and J z . When J z is fixed and D is varied, we carry out 
the FSS analysis based on the following equation: 

0: tr (N,D,J^) = N-^dDc-D)^/"*), (13) 

where a = x or z, jf x is the fixed J Z: and D c is the criti- 
cal point. The same equation concerning the string- type 
order parameter was used in Ref. [l7T ] . The exponents r] a 
and v a arc defined as 

£(N ■ = oo,D,J**)~(D c -D)- v , (14a) 
0^ T (N,D c ,^)^N-^, (14b) 

where £(N,D,J Z ) represents the correlation length. 
When D is fixed and J z is varied, on the other hand, 
the FSS formula is given by 

0£ r (iV,£> fix , J z ) = N-^((J C Z - J^N 1 /""), (15) 

where the critical exponents r\ a and v a are given by 

£(N = oo, D R *, J z ) ~ ( Jl - J)-", (16a) 
0^(N : D R \J c z )~N-^. (16b) 

Note that as the BKT transition point is approached, 
only Eq. pTb|) or Eq. (fl6b]> is realized. In this case, 
the correlation length grows exponentially ; the depen- 
dence is different from Eq. (|14ap or Eq. (|16a|) , having a 
finite exponent. 



Successfully obtaining a universal function "J irrespec- 
tive of the system size N in the critical region near the 
critical point, allows us to determine the critical point 
and critical exponents. However, the FSS analysis still 
has a problem in that the width of the critical region 
is unknown. Since the width depends on the values of 
the control parameters, it is difficult to determine or es- 
timate an appropriate width. To avoid this difficulty, we 
perform the extrapolation explained below. 

C. Ground-state 
phenomenological-renormalization-group analysis 

In this paper, we perform a procedure to obtain 
the transition point and the critical exponents con- 
sistently considering the ground-state quantities. We 
call this procedure the ground-state phenomenological- 
renormalization-group (GSPRG) analysis. The first step 
of the procedure is to examine a finite-size quantity de- 
fined as 

n M log [Q^(N 2 ,D, J Z )/Q^(N 1: D, J z )} 

rj a {N,D,J z ) = — , 

log[7Vi/7V 2 J 

(17) 

where N = (Ni + N 2 )/2, N 2 = AN + N u and AN = 6. 
Here, we examine the direction a such that Of tl . shows 
critical behavior. The quantity indicates the gradient of 
the curve of the dependence of 0" tr on 1/N in a plot 
with both the axes on the logarithmic scale. The gradi- 
ent should be constant for large system sizes when the 
set of D and J z corresponds to the boundary of the 
Haldane phase. The quantity rj a (N, D, J z ) converges to 
the critical exponent r\ a defined by Eq. (|14b[) or Eq. 
(|16b|) for (D, J z ) on the boundary; the N dependence 
of r] a (N, D, J z ) shows a stable convergence to a finite 
value when the system size is increased. On the other 
hand, when the point (D, J z ) is not on the boundary, 
r] a (N, D, J z ) shows a different behavior. For (D, J z ) in- 
side the Haldane phase, 0" tI tends towards a nonzero 
value as N is increased. Thus the gradient rj a (N,D, J z ) 
rapidly decreases. For (D, J z ) outside the Haldane phase, 
0" tr decays rapidly with increasing N. This decay is more 
rapid than that for Eq. (|14b[) or Eq. (|16b|) . In this case, 
the gradient rj a (N,D,J z ) rapidly increases. Therefore, 
we can find the critical point from the difference in the TV- 
dependence of rj a (N,D, J z ). The difference is expected 
to be more apparent when N increases sufficiently to di- 
minish the edge effect. In order to estimate the critical 
point, we have investigated the behavior of rj a (N, D, J z ) 
for N = 6, 12, • • • , 90, 96 and found the characteristic be- 
havior of rj a (N, D, J z ) in the region of large N. The 
numerical procedure to determine the critical point by 
observing the behavior of rj a (N, D, J z ) is as follows. 

• When the differentiation of the finite-size quantity 
r] a (N, D, J z ) satisfies the following conditions for a 
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large system size, 

r]W(N,D,J z ) >OA V W(N,D,J z )>0, 



V W(N,D,J Z ) <OAT)W(N,D,J z )<0, 



(2) 



(18a) 



(18b) 



we can consider that the system with (D,J Z ) is at a 

critical point. Here rj^ (N , D, J z )(n — 1,2) is the 
numerical differentiation given by 

d 



r,W(N,D,J z ) = 



d(l/N) 



Va (N,D,J z ). (19) 



The differentiation is approximated by the differ- 
ence because N is integer or half integer. 



If the differentiation reveals 

77« (N, D, J z ) > A r)W(N, D, J z ) < 0, 



(20) 



we can consider that rj a (N,D,J z ) will decrease 
rapidly with increasing system size. In this case, 
the system with (D, J z ) is in the ordered phase with 
respect to the string order. 



If the differentiation satisfies 

r/« {N, D, J z ) < A r,^(N, D, J z ) > 0, 



(21) 



we can consider that ri a (N,D,J z ) will increase 
rapidly with increasing system size. In this case, 
the system with (D, J z ) is in the disordered phase 
with respect to the string order. 

We summarize the difference in the behavior of 
r] a (N,D, J z ) in Table [IJ Now, the boundary for a finite- 
size system between the critical region and the string- 
ordered phase is given by N' defined in 



7? W (jv' > D , j z ) > o A V M (N' , D , J z ) = 0, 
or 

r)£\N', D, J z ) = A 7${N', D, J z ) < 0, 



(22a) 



to find a boundary between Eq. (|18a[ I18b[) and Eq. 
([2H)l . We obtain N' as a real positive number because 
r]a\N', D, J z ) is an interpolated value of r]a\N,D,J z ) 
for i = 1,2. The boundary between the critical region 
and the string disordered phase is, on the other hand, 
given by N' defined in 



n£ ) (N',D,J z ) = QAr ] W(N',D,J z ) > 0, 
or 

77 W (N\ D, J z ) < A i$ (JV', D, J z ) = 0, 



(23a) 



(23b) 



to find a boundary between Eq. (|18a[ Il8b[) and Eq. 
pTjl . When the critical behavior appears only at a 
point, the width between the critical-ordered boundary 
and the critical-disordered boundary shrinks as N' in- 
creases. Such behavior will be presented in section HV Bl 



TABLE I: Behavior of r/ a (N , D, J z ) in the four phases and 
near the three transition lines. HN, HL, HX: Haldane-Neel, 
Haldane-Large-D, Haldane-XY. RD, RI, SF: rapidly decreas- 
ing finite-size exponent (JT7J, rapidly increasing exponent, sta- 
bly finite exponent, as the system size is increased. 





Haldane 


HN transition line 


Neel 


T] X (N,D,J Z ) 


RD 


SF 


RI 


r, z (N,D,J z ) 


RD 


RD 


RD 




Haldane 


HL transition line 


Large-D 


Va (N,D,J z ) 


RD 


SF 


RI 




Haldane 


HX transition line 


XY 


Va (N,D,J z ) 


RD 


SF 


SF 



In this work, we consider the width of the critical region 
between the two boundaries to be an error in our analysis 
of the transition point if the width is very narrow. In a 
case of the BKT transition, the critical-disordered bound- 
ary does not appear. In this case, we must estimate the 
transition point carefully only from the critical-ordered 
boundary. Details of this treatment will be given in sec- 
tion [EvEl 

The estimation of r/ a of the string order parameter by 
PRG analysis has been reported by Hida [2(| based on 
finite-size data of the string orders by an exact diagonal- 
ization (ED) method. Since the system size is limited to 
being very small, however, the finite-size effect becomes 
significant. To avoid this difficulty as much as possible, 
Hida combined the critical point determined from the 
energy gap under the open boundary condition and the 
string correlation functions under the periodic bound- 
ary condition. For our purposes, we impose only the 
open boundary condition for our DMRG calculations and 
employ a definition of the string order Eq. (TTj) as the 
longest-ranged component. 

In the final stage of the present analysis, we estimate 
the critical exponents v a (D c , J^ x ) and ^ Q (D fix , J^). We 
first consider the case where D is controlled for a fixed J z . 
Within the FSS analysis based on Eq. (fTB")) . an appropri- 
(22b) ate set of D c , rj a , and v a is expected to give a universal 



function ^ near D — D c independent of N. However, it 
is difficult to determine the width of the critical region, 
as we have mentioned. We instead focus our attention on 
the gradient of the universal function * at D = D c . We 
note that the iV-independence of the gradient is a neces- 
sary condition for the existence of the universal function 
near D = D c . Therefore, we assume that the gradient 
for N = Ni and that for N = N2 agree with each other 
for the same D c , rj a , and v a , to give 



dD 

d[OUN2,D,J 



D=Dc 
fixM 



dD 



(24) 



D=Dc 



We input D c and r] a (N, D, J z ) determined above into D 
and r\ in this equation and solve with respect to \jv a . 
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Denoting the solution by l/v a (N, D c , Jf x ), we obtain 
1 



r, a (N,D c ,J^ 



Iix\ 



v a (N,D c ,J**) 

log[(OQ s " tr (7V 2 , D c , J**)/(0>)% tI (N 1 ,D c , J^)) 
\og[N 2 /N x ] 

where (0')" tI (N, D c , J z x ) represents 

fix , dO« tI (N,D,J^) 



(OXMD c ,J: 



3D 



(25) 



(26) 



D=D C 



We note that v a (N, D c , J^ x ) is a finite-size quantity and 
we examine the A^-dependence of this quantity. An ex- 
trapolation to the limit N — > oo provides the expo- 
nent v a (D c , Jf x ). Hereafter, we call rj a (N,D,J z ) and 
v a (N, D c , J^ x ) the finite-size exponents. 

We next consider the case where J z is controlled for a 
fixed D. The same derivation as the above from Eq. (fT5|) 
leads to the finite-size exponent: 



1 



77 a (7V,D fix ,J=) 



\og[{Q')^{N2,D^,J c z )/{OXtANuD^,^) 
\og[N 2 /N x ) 

where (0')" tr (N, D fix , Jf) represents 

dO? tr (N,D R *,J z ) 



(0% I (N,D»*,J c z ) = 



(27) 



(28) 



The extrapolation of v a {N , Z? fix , J z ) gives the exponent 
v a {D^,Jl). 
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FIG. 1: (D, J z ) dependence of the order parameters for the 
Hamiltonian (TT} on (a) the J z = 0.5 line and (b) the D — 
line at N = 300. +: Oj w , x: Of tr , □: C^ el , Q= Ofc. 



IV. RESULTS AND DISCUSSIONS 
A. Behavior of order parameters 

Let us first review the behavior of the four order pa- 
rameters under consideration in a finite-size system and 
summarize some important relations between them. In 
a moderately large system, we can see indications of 
asymptotic behavior in each order parameter, although 
slow convergence prevents a full description. Some are 
characteristic for a given region of the parameter space, 
which is specified as one of the Haldane, Neel, Large-D, 
and XY phases. 

We illustrate O b a tr (300, D, J z ) and O^ 6cl (300, D, J z ) 
with a = x or z in Fig. [T] The D-dependences of 
the order parameters on the line J z — 0.5 and their J z - 
dependences on the line D = are shown in Fig. [IJa) 
and QJb) , respectively. We now compare the magnitudes 
of the four order parameters, (i) If J z > 0, 0" tI is larger 
than O^eci- This is a known relation found by Kennedy 
and Tasaki [2l|. (ii) When D = and J z = 1, we have 
Og tl . = Og tr due to the isotropy of the system, (iii) When 



J z is decreased and when J z crosses a critical point at 
J z = 0, Og tr is smaller than 0^ 6cl . This fact will be dis- 
cussed and utilized in §IVDI Although 0^ 6cl appears to 
be nonzero around < D < 1 with J z = 0.5 and around 
< J z < 0.5 with D = 0, we can confirm that 0§ 6el 
in this region vanishes for the long-ranged limit. On the 
other hand, OJ tr around < J z < 0.4 with D = looks 
very small but it survives as a nonzero quantity in an 
infinite system, as shown in i jlVI 

The phase boundaries of Haldane-Neel, Haldanc- 
Large-D, and Haldane- XY are denoted by dotted lines, 
though they are given only as indicators as we will de- 
termine the boundaries in the following subsections. We 
can see that some or all of the order parameters van- 
ish at the phase boundaries. Also, the order parameters 
are continuous around the boundaries, which suggests 
that the phase transitions are continuous. Therefore, the 
FSS analysis and the GSPRG procedure are feasible for 
capturing critical phenomena in this case, except for the 
Berezinskii-Kosterlitz-Thouless (BKT) transition which 
does not satisfy the conditions of Eq. (|I4a[) or (|16a|) and 
thus requires extra consideration. For GSPRG, however, 
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it is possible for us to capture the transition by looking 
at the exponent rj a as discussed in i )IVD[ 

We next observe the behavior of the four order param- 
eters in each phase to determine their thermodynamic 
limits. In Fig. [21 we illustrate the behavior of the order 
parameters (jTTJ) as a function of the inverse of the sys- 
tem size at the representative points (D, J z ) = (—2,0.5), 
(2,0.5), (0,0.5), and (0,-0.5). These sets of parameters 
correspond to the Neel phase, Large-D phase, Haldane 
phase, and XY phase, respectively. We observe that in 
the Neel phase, only the order parameters in the z direc- 
tion remain nonzero in the limit N — > oo. All the four 
parameters vanish in the Large-Z) phase in the thermo- 
dynamic limit. We note that in the Haldane phase, only 
the string order parameters in the two directions remain 
nonzero in the thermodynamic limit. It is difficult to 
judge in Fig. [DJd) whether both of the transverse order 
parameters in the XY phase vanish or remain nonzero 
in the thermodynamic limit. We plot the same data on a 
logarithmic scale in Fig. [3J For large N, the data exhibits 
a linear behavior, which suggests that the transverse or- 
der parameters in the XY phase are critical, consistent 
with previous reports [1, [l6| . Consequently we can con- 
firm that the order parameter O q (oo, 0,-0.5) vanishes in 
the XY phase, in the limit N — > oo. 



B. Haldane— Large-D transition line 

In this subsection, we examine the transition from the 
Haldane phase to the large-!? phase. This transition is 
known to be of Gaussian type. As we observe in ^IV Al 
the string order parameters in the Haldane phase remain 
nonzero for both a — x and z while both of the Neel order 
parameters vanish along the directions x and z. We also 
observe critical behavior near the transition line in both 
O'itr and O s V 

To begin with, we consider difficulties in the FSS anal- 
ysis near the transition between the two phases. In this 
analysis, we have adjusted the critical point D c and ex- 
ponents rj a and v a so that the data for N = 24, 48, and 
96 follows a universal function. The results are depicted 
in Fig. [U In Fig. Ufa), we observe a deviation from the 
universal function ^ at D, not far from D = D c . The 
appearance of this deviation depends on the system size 
and the direction a. Thus, it is not easy to determine 
the critical region around D = D c with finite-size data 
less than 100 sites in this case. 

Despite this difficulty, we can choose input parameters 
D c , i/ a , and r) a such that a universal function ^> appears 
near the transition point. In Fig. Ufa), the string corre- 
lation functions in the direction a — x provide us with 
£>H L (J Z = 1.25) = 1.16 and v x {D = 1.16, J, = 1.25) = 
1.05. On the other hand, in Fig. 2{b), the string cor- 
relation functions in a = z give D^ h (J z = 1.25) = 1.16 
and u z {D = 1.16, J z = 1.25) = 1.20. The estimate of the 
transition point D c for a — x and that for a = z agree 
with each other. This fact strongly suggests that the 
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FIG. 2: Order parameters for the Hamiltonian (flj as a func- 
tion of the inverse of the system size. Each panel represents 
(a) the Neel phase, (b) the Large-D phase, (c) the Haldane 
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FIG. 3: Order parameters for the Hamiltonian |T} in the 
XY phase as a function of 1/N. A fitting based on O(N) = 
dN' 1 is carried out for N = 96, 90, • • • , 72. □: 0% M , O 

^str ■ 



string correlation functions for the transverse and lon- 
gitudinal directions reveal a common phase transition. 
We should note that 0^ tI and 0| tr are clearly different 
quantities near the transition point, because there are dif- 
ferences in their exponents, for example r\ x — 0.312 and 
?y 2 = 0.756. Our FSS analysis gives u x (D — 1.16, J z = 
1.25) = 1.05 and v z (D = 1.16, J z = 1.25) = 1.20. We 
recall that the growth of the correlation length deter- 
mines the critical behavior near the transition point from 
a general argument of the renormalization group concern- 
ing critical phenomena. In this framework, only a single 
characteristic length in a system shows critical behav- 
ior. The characteristic length must be the correlation 
length of the system. Thus, the exponent of the corre- 
lation length should be unique for the order parameters. 
In this case, the correlation functions of the string order 
parameters along both a = x and a = z show critical 
behavior as shown by the FSS analysis. From this argu- 
ment, v x and v z should exhibit a serious finite-size effect, 
which we will examine and solve by GSPRG analysis. 

We consider the case of J z = 1.25 in order to ob- 
serve the finite-size effect. In Fig. [5] we illustrate our 
results for the exponents r] x (N,D, 1.25), r] z (N,D, 1.25), 
v x (N, 1.16, 1.25), and u z (N, 1.16, 1.25) determined by 
GSPRG analysis. In Figf5ja), we observe the critical- 
disordered boundary at D — 1.17,7V' = 55.3. The 
critical-ordered boundary is also observed at D = 
1.15, N' = 45.5. On the other hand, we obtain no bound- 
aries defined in eqs. (j2^ - (j23bf in the case of D = 1.16. 
This fact suggests that the critical region for finite-size 
systems in our study is realized around D = 1.16 with 
a narrow width. In order to confirm whether the width 
shrinks or not as the system sizes increase, we examine 
the relationship between N' and D so that the case is on 
the boundary. We obtain some of the critical-disordered 
boundaries at (D = 1.19, N' = 36.8), (D = 1.2, N' = 
32.7), and (D = 1.21, N' = 29.8). We also obtain some of 
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FIG. 4: FSS of (a) the transverse string order and (b) the 
longitudinal string order on J z — 1.25 near the Haldane- 
Large-D critical point. Q, +, and x represent N = 24, 48, 96. 



the critical-ordered boundaries at (D = 1.14, N' = 32.7), 
[D = 1.13, N' = 30.0), and (D = 1.12, N' = 29.1). These 
results indicate that the critical region for a given N' is 
gradually narrower when N' increases although the ex- 
pression of the relationship between N' and D on the 
boundary is unknown in the present stage. It is reason- 
able to conclude that the critical region shrinks and goes 
to the transition point for the infinite-size system. When 
one can confirm whether the critical region between the 
two boundaries is sufficiently narrow or not, the width of 
the region should be regarded as an error coming from 
the maximum system size and the interval of D in the 
performed calculations. In this work, thus, we conclude 
Df L {J z = 1.25) = 1.16 ± 0.01. Note here that we can 
obtain the same critical point from r\ z in Fig[5jb) in the 
same manner. Hereafter, we determine critical points 
with an error in this way. In order to confirm whether the 
critical behavior (|14b[) or (|16b[) appears or not in the orig- 
inal correlation functions, each string correlation function 
as a function of 1/N is shown in the logarithmic scale in 
inset figures. The finite-size string correlations for each 
direction clearly reveal a power-law decay behavior at the 
critical point D^ L (J Z = 1.25) = 1.16. On the other hand, 
a behavior deviating from power-law decay appears in 
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FIG. 5: GSPRG analysis of (a) r} x , (b) r) z , and (c) u a (D = 
1.16) on J z = 1.25. N is given by Eq. (JTZJ. The dotted 
curves in (a) and (b) are guides for the eyes. The dotted lines 
in (c) are linear fitting lines. The string correlation functions 
0* tI {N,D, J z = 1.25) and 0% tt {N,D, J z = 1.25) as a function 
of 1/N are shown in the logarithmic scale in the inset figures 
of (a) and (b), respectively. 



the cases of D = 1.0 and D = 1.3 in the ordered and 
disordered phases, respectively, as we have mentioned 
in SjTTTCj Note here that a comparison with these in- 
sets shows that the system size dependence of the finite- 
size quantity (jTTJ) sensitively change near the transition 
point. We next observe the N dependence of the finite- 
size exponents of v x (N , 1.16, 1.25) and v z (N, 1.16, 1.25) 



for J z = 1.25 and Df h (J z = 1.25) = 1.16 in Fig^c). 
These two finite-size exponents, v x and v Zl get gradu- 
ally closer with increasing N . In the limit N — > oo, v x 
and v z appear to approach a single value ^1.2. This is 
consistent with the above argument on the unique char- 
acteristic length. Consequently, the problem of the dis- 
agreement of v x and v z in the FSS analysis occurs due to 
the finite-size effect and is resolved by GSPRG analysis. 

We now consider the transition point D c for a fixed 
J z = 1. In this case, many studies have reported various 
estimates for the boundary of the Haldane phase, D c : 
D c = 0.93 ± 0.02 in Ref. U^JD C = 0.99 ± 0.02 in Ref. 
[23, D, = 0.90 ± 0.05 in Ref. J0, D c = 1.001 ± 0.001 in 
Ref. H, D c = 0.95 in Ref. (24^ D c = 0.99 in Ref. [3, 
and D c ~ 0.97 in Ref. [26|]. Among these works, only a 
single study [15j was based on the analysis of the string 
order, although data from the numerical-diagonalization 
calculations in this study for small clusters might not be 
sufficient to show the transition point. Recently, Tzeng 
and Yang [2fJ| investigated the fidelity susceptibility (23| 
of the ground state by the DMRG method to detect quan- 
tum phase transitions for the system. This work exam- 
ines only the information of the ground state, a feature 
that is shared with our present analysis. Other works 
analyzed the structure of low-energy levels. From the 
present analysis, our estimate is D c = 0.975 ± 0.015, 
which we have obtained irrespective of a = x or a = z. 
Although the estimates are all very close to each other, 
there are small differences between them even taking er- 
rors into account. The reason for these differences is not 
clear at present and should be resolved as a future issue. 

Next, we consider the transition point D c for a fixed 
J z = 0.5. The estimation of this point is suitable for 
checking the availability of our analysis procedure, be- 
cause a relatively large exponent v which is reported 
2.38 by analysis of the energy level structure appears 
[l(J. Several previous studies presented numerical data 
of the transition point as follows: D c = 0.635 in Ref. 
Ik A: = 0.65 in Ref. [l(j, D c = 0.633 ± 0.02 in Ref. 
28|| . All of these works examined the free energy near 
the critical point to determine the critical point. In par- 
ticular, the recent study [H| develops rapidly converging 
methods by using the differentiations of a quantity, which 
is derivative of the ground state energy with respect to 
a controlled parameter, as a function of TV. From the 
viewpoint of using only information in the ground state 
for detecting a quantum phase transition, our analysis 
and their analysis have a common policy. Our estimate 
for D C (J Z = 0.5) is 0.67 ± 0.04, and this estimate is also 
consistent with all previous reports. 

In accordance with the above results, we apply the pro- 
cedure to estimate the critical behavior for other J z , con- 
firming the J z dependence of D^ L (J Z ) and v(D^ L , J z ). 
The error of v(Df Jj ,J z ) is estimated by \v{D^ h ,J z ) - 
V ( D N =93,HL ; We illustrate our results in Fig. H to- 
gether with those of previous reports [H, [ll], [25[. Our 
estimates of D^ L (J Z ) and v(D^ L , J z ) are common for 
a — x and a — z within errors. Our transition line is al- 
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d^J), this value of Vn6 C i(= 1 / 7 ?Neei) S ives v ~ L33 > which 
is consistent with our estimate shown in the inset of Fig. 
|B1 This consistency also supports the scaling hypothe- 
sis that the growth of the unique correlation length de- 
termines all the critical behavior around the transition 
point. 




FIG. 7: 1/N dependence of 0^ kel {N,D, J x ) at J z = 1 and 
D c = 0.975. The dotted line shows 0^ 6cl {N, D, J z ) oc N~ v ^ 

With T^eel = 0.40. 



FIG. 6: (a) Haldane-Large-D transition line and (b) critical 
exponent v a {Jz)- : evaluated value of x-component, 0: 
evaluated value of z-component, +: Ref. □: Ref. [llj], A: 
Ref. (25|. The inset of (b) magnifies the data at J z — 1 to 
allow a clear comparison to distinguish the data. 

most consistent with those of previous reports [t| [ll|, [25| , 
in which the energy-level structure is analyzed. Our es- 
timates of v also agree well with previous reports within 
errors. Consequently our GSPRG analysis successfully 
captures the transition between the Haldane phase and 
the large- D phase. 

The correlation length exponent v is known to be re- 
lated to other critical exponents. In the Gaussian tran- 
sition, Okamoto obtained the following relationship from 
the argument by the bosonization method: 

V =T~^> (29) 

where ??N6ci 1S exponent defined by (SqS^) ~ 
(— l) r r~' ? N«ei at the transition point. Note here that 
^Neci^Neci = 1 holds. To confirm the consistency be- 
tween our estimate of v and the decay of the Neel corre- 
lation function, we plot our 0^ €el (N, D, J z ) at J z = 1 
and D c = 0.975 as a function of 1/N on a logarith- 
mic scale in Fig. [7] We clearly observe a linear be- 
havior for large N. We have added the dotted line 



C. Haldane— Neel transition line 

In this subsection, we examine the transition from the 
Haldane phase to the Neel phase. This transition is con- 
sidered to be of Ising type. We recall that in a transition 
of Ising type, the exponent of the correlation length is 
v = 1 when the system approaches the transition point. 

We have mentioned in the above that the longitudinal 
string order is nonzero in both of the Haldane phase and 
the Neel phase and that the order does not reveal the crit- 
ical behavior at the transition point. This means that the 
longitudinal string order is not appropriate for studying 
the Haldane-Neel transition. Therefore, to study this 
transition we examine only the transverse string order. 
By GSPRG analysis of this order, we determine the tran- 
sition point ,J C Z for a given D or the transition point D c for 
a given J z and the critical exponent v near the transition 
point. 

We consider the case of D = 0.5. We illustrate our 
result for finite-size exponents r/ x and v in Fig. [S^a) 
and (b), respectively. Our estimates are J Z ' W (D — 
0.5) = 1.4905 ± 0.0015 and v( J z = 1.4905, D = 0.5) = 
1.006 ±0.016. Our estimate of the transition point is dif- 
ferent from that of J Z ' HN (D — 0.5) = 1.536 reported in 
Ref. 0] . To find the reason for the difference between the 
two estimates, we have made numerical-diagonalization 
calculations of finite-size clusters up to N = 22 under the 
periodic boundary condition and obtained the eigenener- 
gies of the low-energy states. We have performed the 
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for D = 0.5 from our DMRG results for AiV = 6. The dotted 
curves in (a) are guides for the eyes. The dotted line in (b) 
is a linear fitting line applied to the data for large N . (c) 
Extrapolation procedure of the finite-size critical point Jl for 
D — 0.5 by the numerical-diagonalization method for AiV = 
2. The numerical-diagonalization data for N = 17, 19, 21 are 
new in this paper. The dotted line is the extrapolation line 
in Ref. 0. The broken curve is a guide for the eyes. 



same analysis as that in Ref. [9j and determined the 
finite-size critical point J^ ,HN (iV, D — 0.5) as J z at which 
the scaled energy gap does not depend on the system size 
for AiV = 2. The results are depicted in Fig. [H[c) . From 
our numerical data for N — 8, 10, • • ■ , 16, we successfully 
reproduce the results of Ref. Q. On the other hand, 



we can observe that JJ< HN (N, D = 0.5) of N = 17, 19, 21 
gradually departs from the fitting line of the extrapo- 
lation in Ref. @. Our new data points approach our 
estimates from the string order by the DMRG calcula- 
tions, as shown by the guide for the eyes denoted by 
the broken curve in Fig. |UJc). This agreement suggests 
that the results from the numerical-diagonalization and 
DMRG calculations are consistent with each other if we 
accept the interpretation suggested by the broken curve. 
Hence, careful extrapolation with respect to system size 
is required. 

The N dependence of our new data appears exponen- 
tial rather than polynomial. A similar N dependence of 
jc.HN^y-^ D = 1) was reported in Ref. [TlJ] , in which cal- 
culations up to N = 48 based on the multi-target DMRG 
method with an infinite-system algorithm were carried 
out under the periodic boundary condition. Our result 
and Ref. [ll| suggest that the absence of polynomial 
components does not depend on the values of the param- 
eters of the system. It is important to be careful when 
a system-size extrapolation of an Ising transition point 
is carried out by the PRG analysis of the energy-level 
structure. 

We now compare estimates of the transition point be- 
tween Ref. [llj and the present analysis. Reference 11 1 
gives J^' HN (Z3 = 0) = 1.186. From the present analy- 
sis of our data up to N = 150, we obtain J c z ' HN (£) = 
0) = 1.1860 ± 0.0003 for the transition point. Our esti- 
mate, with a very small error, agrees excellently with the 
estimate in Ref. 
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We now discuss our estimate of v. Our estimate v(J z = 
1.4905, D = 0.5) = 1.006 ± 0.016 is in good agreement 
with v — 1 of the Ising-type transition. This agreement 
also suggests that our analysis successfully captures the 
Haldane-Neel transition as well as the Haldane-Large-D 
transition. 

We can now summarize our results for the transition 
points D c for a given J z and the critical exponents v be- 
tween the Haldane and the Neel phases from our DMRG 
data. The results are depicted in Fig. [SJ Figure [5Ja) 
shows that our estimates for the transition points are 
in good agreement with the results in Ref. [11| of the 
multi-target DMRG method and the results in Ref. 
of the numerical diagonalizations. In Fig. [DJb), our es- 
timates for the exponent agree well with v = 1 irrespec- 
tive of J z . Note here that the center values of our esti- 
mates, namely the extrapolated results, are much closer 
to v = 1 than the results in Ref. [llfl . although our 
errors are estimated to be larger. Note also that the 
error in v(D™ = 2.015, J z = 2.5) is quite large. The 
reason for this is considered to be that the curve of the 
Haldane-Neel transition points and that of the Haldane- 
Large-Z) transition points approach each other. A sim- 
ilar phenomenon appears when the central charge c on 
the curve of the Haldane-Large-D transition points was 
estimated in Ref. [9j, in which the estimate of c gradu- 
ally deviates from c = 1 around J z > 1. In the report 
of Tzeng and Yang (26|, the transition point and the 
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FIG. 9: (a) Haldane-Neel transition points. : Our work, 
+ : Ref. [jj], 0: Ref. [ill ], (b) Haldane-Neel critical exponent 
v. The inset figure compares our v with that of Ref. [ill ]. 
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OnccI = 0. 



under the condition 
O 

This means that the region 

Ostrl — I^Nccllj 



(31) 



(32) 



is not in the Haldane phase because the inequality 
(I3T))) and Eq. (pTTj) cannot both be satisfied at the 
same time assuming the inequality (|32p . However, 
it is not as easy to make a direct comparison of 
these quantities in the limit N — > 00 as for the in- 
equality (|3"2"|) . We can instead compare the finite-size 
quantities \0£ T (N, D, J z )\ and \0^ 6ol (N, D, J z )\. Re- 
call that for N = 300, \0" tT (N, D, J z )\ is smaller than 
\0^ el (N,D,J z )\ when J z < 0, whereas \0% x (N,D,J z )\ 
is larger than \0^ 6el (N, D, J z )\ when J z > 0. We have 
studied the system size dependence of this behavior; our 
results are depicted in Fig. 1101 The behavior is clearly in- 
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critical exponent are given as Dc HN (J z = 1) ~ —0.31, 
i/(D c HN = -0.31, J z = 1.0) - 1.05, respectively, from fi- 
delity susceptibility analysis. Our estimated values at the 
same point are Dc = -0.315 ± 0.003, v = 1.004 ± 0.019, 
which are more precise than the values of Tzcng and 
Yang. 



D. Haldane— XY transition line 

In this subsection, we examine the transition from the 
Haldane phase to the XY phase. This transition is con- 
sidered to be a BKT-type transition. We recall that in a 
BKT-type transition, the exponents rj x = 1/4 and ri z = 1 
appear at the transition point and the exponent v cannot 
be defined because the correlation length grows exponen- 
tially. 

We consider the case J z < and examine the magni- 
tudes of the string order and the Neel order. We refer 
back to the behavior of orders characterizing the Haldane 
phase, in which we have 

|O s z tr | >0and |0* tr | >0, (30) 



FIG. 10: Order parameters crossing at D = — 1.0. O an d x 
represent 0^ iel and Of tr . 

dependent of system size. We can also confirm this inde- 
pendence irrespective of D for cases between the Haldane 
phase and the XY phase. Our present results suggest the 
inequality (|3"2"|) and indicate that the Haldane-XY transi- 
tion point satisfies J z ,iiXY (D) > 0. The finding is entirely 
consistent with previous works. Thus, it is sufficient to 
consider the case of J z > hereafter in examining the 
Haldane-XY transition. 

We now estimate J Z ' HXY by our GSPRG analysis. We 
consider the case D — —0.5 for J z > 0. For this pur- 
pose, we examine the finite-size exponent r) a (N,D — 
— 0.5, J z ), and estimate the critical-ordered boundary 
point ri { a\N',D = -0.5, J z ) given by Eq. |23a|) or Eq. 
(|23b[) . Our results are depicted in Fig. [TTJ We find that 
the critical-ordered boundary given by Eq. (|23b[) appears 
when J z is 0.1, 0.15, 0.18, and 0.2, but it does not appear 
when J z is 0. Concerning with x-component of the string 
order, we find the boundaries at ( J z = 0.18, N' — 29.6) x , 
(J z = 0.15, N' = 31.0)3., and (J z = 0.1, N' = 34.6)*. 
Concerning with z-component of the string order, on the 
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FIG. 11: GSPRG analysis of the exponent r] a {N,D = 
—0.5, J z = 0). The dotted lines are the rigorous exponent 
values of the BKT transition: r\ x — 0.25, rj z = 1. Oj +j x,D, 
and represent J z = 0, 0.1, 0.15, 0.18, and 0.2, respectively. 



other hand, we have (J 2 = 0.18, N' = 32.6) z , (J z = 
0.15, N' = 35.2),, and [J z = 0.1, N' = 41.3), as the 
boundaries. In the cases of both of the components, one 
can observe that N' grows when J z approaches J z = 0. 
These phenomena lead to our result that J^ ,HXY is be- 
tween J z — and J z — 0.1. For estimating the transi- 
tion point more accurately, the critical-ordered boundary 
point is extrapolated to the limit N' — > oo. The results 
are depicted in Fig. [T2] Since the leading dependence of 
jc, hxy on i/jy> j s unknown, we here choose the power 
1/N' so that the dependence is almost linear. We can 
successfully determine an appropriate value of the power 
for each a = x and a = z, although the a — x and a — z 
values differ from each other. A linear extrapolation 
gives J 2 C ' HXY (L» = -0.5) = 0.00 ± 0.10 from the trans- 
verse component and J^ HXY (L> = -0.5) = 0.01 ± 0.08 
from the longitudinal component. Here we determine 
the error as being the difference between the values ob- 
tained by the extrapolation and the finite-size critical 
point J^' HXY for maximum N'. Both results suggest 
jc,hxy^£) _ _o.5) ro irrespective of the direction of 
the string order parameter, which is consistent with a 
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FIG. 12: Behavior of the critical-ordered boundary point 
jc, hxy a j. j-y _ _q g as a f unc tion of the system size N' . 
Results are shown for tuf = 150. The error is estimated 
from the difference between the results for tuf ~ 100 and for 
rriF = 150. 



previous report Q . 

Next, we examine what type of transition this is. Our 
finite-size exponents in Fig. [Tl]at J z = indicate rj x (D = 
-0.5, J z = 0) ~ 0.25 and rj z (D = -0.5, J* = 0) ~ 1.0. 
These values agree well with the exponents of the BKT 
transition r/ x = 1/4 and rj z = 1. Our results are also 
consistent with many previous works [1, d, H, [3(| • There- 
fore, our GSPRG analysis applied to the string correla- 
tion functions is useful in capturing BKT transitions. 



V. SUMMARY AND REMARKS 

We have investigated critical behavior near the bound- 
ary of the Haldane phase in the ground state of an 
anisotropic S — 1 chain from the viewpoint of string 
correlation functions estimated precisely by standard 
finite-size DMRG under the open boundary condition. 
We have developed the ground-state phenomenological- 
renormalization-group analysis and used it to analyze the 
correlation functions. This analysis provides us with the 
transition point of the boundary of the Haldane phase 
and the critical exponents at and near the transition 
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point. Our estimates for these quantities agree with 
those previously obtained from analysis of the energy- 
level structure. 
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FIG. 13: Phase diagram for the Hamiltonian (TTJ. The 
Haldane-Large-D, Haldane-Neel, and Haldane-XF transi- 
tion points are denoted by 0> A, 0, respectively. All the lines 
are guides for the eyes. x-Haldane and z-Haldane represent 
the Og tr > O z str region and Of tf < Of tr regions, respectively. 

We summarize the transition points as a ground-state 
phase diagram in Fig. [T3j This figure presents the 
phase boundary of the Haldane-Large-Z?, Haldane-Neel, 
and Haldane- ATU transitions. Note additionally that the 
dominant order parameter is 0^ tI in most of the Haldane 
phase. 

A feature of our approach, GSPRG analysis, is that 
only common quantities under the same condition are 



treated in a unified manner irrespective of the type of 
phase transition. Although we have employed the DMRG 
method in this paper to calculate the order parameters, 
we are not limited to the DMRG method if we can obtain 
precise estimates of the order parameters. The string 
order parameters of the Haldane phase in the S = 1 
chain are examples. Other multi-point correlation func- 
tions for finite-size clusters may also be applicable. If we 
precisely calculate an appropriate ground-state quantity 
that plays the role of an order parameter, the framework 
of the analysis would be widely applicable for capturing 
ground-state critical behavior irrespective of the method 
of calculation and the kind of order parameter. We hope 
that the procedure presented in this paper contributes to 
future studies of quantum phase transitions. 
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